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A Lattice Study of Semi-leptonic Decays of Z)-Mesons 

UKQCD Collaboration - presented by J. Nieves^* 

""Physics Department, The University, Southampton S017 IBJ, UK 

We present results of a lattice computation of the matrix elements of the vector and axial- vector currents which 
are relevant for the semi-leptonic decays D ^ K and D — )• K* . The computations are performed in the quenched 
approximation to lattice QCD on a 24^ x 48 lattice at /3 = 6.2 using an 0(a)-improved fermionic action. 
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1. Introduction and Simulation Details 

There is increasing evidence that quantitative 
calculations of weak decay amplitudes can be ob- 
tained by lattice QCD simulations. D decays pro- 
vide a good test of the method, since the relevant 
CKM matrix elements (Ks and Vcd) are well con- 
strained in the Standard Model. In addition com- 
parisons bctwoon D and D decays reveal the size 
of non-leading terms in the Heavy Quark Effec- 
tive Theory (HQET). The main advantage of the 
lattice technique is that it is based on first prin- 
ciples only and it does not contain free param- 
eters besides the quark masses and the value of 
the lattice spacing, both of which are fixed by 
hadron spectroscopy. Moreover, statistical and 
systematic errors in lattice simulations can be 
systematically reduced with increasing computer 
resources. In this study we present the results 
of a lattice calculation of the matrix elements 
{K*\J>'\D) and {K\J>'\D) (.7'' = .S7''(l - 7,5)c is 
the relevant quark weak current) which contain 
the non-perturbative QCD corrections to the de- 
cays D — > Kl^vi and D — > K*l~^vi. We have used 
an improved quark action, proposed by Sheik- 
holeslami and Wohlert [1], which together with 
an appropriate rotation of the quark fields, leads 
to matrix elements free of discretisation errors 
of order aa" log" (a) . When using this action 
the leading discretisation errors are of the order 
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0{asa) and O(a^) [2]. This "improvement" is 
particularly important here since we are study- 
ing the propagation of quarks whose bare masses 
are aroimd one third the inverse lattice spacing. 
We have analysed 60 gauge field configurations 
at = 6.2 on a lattice of size 24^ x 48. Prop- 
agators were computed at three "light" values 
of the hoping parameter, k = 0.14144, 0.14226, 
0.14262, and at Kc = 0.129 which corresponds to 
the mass of the charm quark. In order to de- 
termine the matrix elements mentioned above, 
we compute three-point correlation functions in 
which the D meson is placed at t = 24, the mid- 
point of the lattice, and symmetrize the correla- 
tors about that point using Euclidean time re- 
versal. We give momenta \pd\ = 0, j^a""^ to the 
Z)-meson and momentum q is inserted at the cur- 
rent. We use the local, improved vector current. 
Results are extrapolated to the physical values 
of the strange (ks = 0.1419(1)) and light quark 
masses {Kcrit = 0.14315(2)). The "light" quark 
propagators were computed using local sources 
and sinks; the "charm" propagator was computed 
using gauge-invariant "smeared" source and both 
local and smeared sinks. Full details of the cal- 
culations are contained in ref. [3]. 

2. Renormalisation Constants Zy and Za 

In this section we discuss the difficulties in de- 
termining the form factors of semi-leptonic D 
K, K* decays due to the presence of discretisa- 
tion errors. In general the matrix elements of an 
"improved" -lattice current, , are related to the 
corresponding continuum matrix elements by: 
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Oia.a) + Oia"^) + . . . (1) 

where the renormalisation constant Zy is inde- 
pendent of the external states i and /, and is 
computable in perturbation theory. For light 
quark masses the discretisation errors are small 
but for the charmed quark this is no longer 
the case. In previous simulations, using Wilson 
fermions (for which there exist 0{a) discretisa- 
tion errors not present in this simulation), these 
effects were modelled by using effective (mass- 
dependent) renormalisation constants Zy^-^ and 
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However we wish to stress that the dis- 



cretisation errors are in general different for ma- 
trix elements of currents with different Lorentz 
indices, between different states, and they could 
have a different q'^ dependence than the form 
factors themselves. Formally there is no reason 
to believe that the discretisation errors can be 
absorbed into universal effective renormalisation 
constants. In ref. [3] the size of the lattice arti- 
facts in the present study is discussed, and the 
conclusion is that although these errors are sub- 
stantially reduced by the use of the improved ac- 
tion, nevertheless even in this case we believe that 
they lead to uncertainties of the order of 10% in 
the form factors. In spite of the above discussion 
we have also assumed that they can be modelled 
by effective renormalisation constants Zy-^ and 
Z"^-^ (at least part of the errors can be so ab- 
sorbed). 

3. Results 

In Figs. 1 and 2 we show the form factors /+ 
and Ai as a function of for the semi-leptonic 
decays D ^ K, K* . In both cases the soHd 
line corresponds to the comparison of the pairs 
((j-^, form factor) with the (j^-dependence of the 
form factors determined from a two-parameter 
(form factor at g'^ = and pole mass) pole dom- 
inance fit to our data. Crosses correspond to the 
form factors at = (up to factors Zy-^ and 
Z^-^ respectively) determined in this way. In 
the case of /+ we also compare (dashed Hne) the 
(/^-dependence of our data with that determined 
from a one-parameter pole dominance fit, where 
we fix the pole mass to that of the vector-meson 



D* , which is the resonance in this helicity channel 
and which is also in an excellent agreement with 
the experimental determination of the q"^ depen- 
dence of /+(g2) by the CLEO Collaboration ([4]). 
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Figure 1. f^{q^) as a function of q^a^. 
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Figure 2. Ai(g^) as a function of q^a^. 

For the axial form factor, Ai, the dashed line 
also corresponds to a one-parameter pole domi- 
nance fit, where the pole mass is fixed to 2.5 GeV 
corresponding to the mass of the axial-vector res- 
onance {Dsi) predicted in this helicity channel by 
quark models. In both figures diamonds corre- 
spond to the form factor at = (up to factors 



Zy-" and Z'^-' respectively) determined by using 
one-parameter pole dominance fits. The fits using 
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Tabic 1 



Form factors at cp = 


for the semi-lepti 


onic decays D — ?• K and D ^ K 












Ai(0) 


y(o) 


^2(0) 


Exp. 


World Ave. [4] 
World Ave. [5] 


0.77 ±0.04 
0.70 ±0.03 


0.61 ±0.05 


1.16 ±0.16 


0.45 ± 0.09 


13 = 6.2, SW-O(a) 
= 6.4, Wilson 
/3 = 6.0, Wilson 
= 6.0,5.7, Wilson 
(3 = 6.0, Wilson 
/3 = 6.3, Wilson 
/3 = 6.0, SW-O(a) 


This work 
ELC [6] 

LMMS [7] 
BKS [8] 
BG [9] 
WU [10] 

APE [11] 


0.67 tl 
0.60 ti= tl 
0.63 ±0.08 
0.90 til 
0.73 ±0.05 
0.76 ±0.15 
0.78 ±0.08 


0.70 t^l 
0.64 ±0.16 
0.53 ±0.03 

n s-i +14 +28 

U.OO _j4 _28 

0.66 ±0.03 
0.59 ±0.08 
0.67 ±0.11 


1.01 til 
0.86 ±0.24 
0.86 ±0.10 

-, +45 +48 

1.24 ±0.08 
1.05 ±0.33 
1.08 ±0.22 


0.66 tl? 

U.'lU _28 -4 

0.19 ±0.21 

+14 +24 
u.oy _j4 _23 

0.42 ±0.17 

0.56 ±0.40 

0.49 ± 0.34 



the constrained pole masses also lead to accept- 
able /dof. As can be seen from these figures 
both methods of extracting the form factors at 
= agree remarkably well, which gives us con- 
fidence in our procedure. It can be also seen from 
these figures that the g2 dependence of both form 
factors is reasonably well described by the pole 
dominance model. Qualitatively similar conclu- 
sions are obtained for the rest of form factors: 
f°,Ao,V and A2 (for details see ref. [3] where re- 
sults for the form factors relevant to the decays 
D ^ Tr,p are also given). In Table 1 we show our 
results for the form-factors at = for the semi- 
leptonic decays D ^ K, K*. In our quoted errors 
we have included the uncertainty in the renor- 
malisation constants, discussed in the previous 
section, in order to take into account some of the 
residual discretisation errors. In this table we also 
compare our predictions with the most recent ex- 
perimental world average and with previous lat- 
tice results. Discretisation errors are, in principle, 
larger for Wilson than for improved actions and 
part of the discrepancies between different lattice 
results in Table 1 are due to different values used 
in the literature for the effective renormalisation 
constants Zy-^ and Z"^-^ . We have also found the 
g2 dependence of the form factors in a wide region 
around q^ = 0. Thus we can estimate the phase 
space integrals and obtain the total decay rates 
(Table 2). In both cases (Tables 1 and 2) our 
results are in good agreement with experimental 
data. The analysis of the B ^ Tr,p matrix ele- 
ments is in progress. 



Table 2 



Semi-leptonic widths for D ^ K, K* decays. 





T{K) 


T{K*) 


(f^)K. 


Exp. [4] 
Exp. [5] 
This work 


9.0 ±0.5 

7.1 ±0.6 
7.0 ±2.0 


5.1 ±0.5 
4.5 ±0.5 

6.0 tr. 


1.2 ±0.1 
1.1 ± 0.2 
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